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ABSTRACT 

We present 5-GHz MERLIN radio images of the microquasar GRS 1915+105 dur- 
ing two separate outbursts in 2001 March and 2001 July, foUowing the evolution of the 
jet components as they move outwards from the core of the system. Proper motions 
constrain the intrinsic jet speed to be > 0.57c, but the uncertainty in the source dis- 
tance prevents an accurate determination of the jet speed. No deceleration is observed 
in the jet components out to an angular separation of 300 mas. Linear polarisation 
is observed in the approaching jet component, with a gradual rotation in position an- 
gle and a decreasing fractional polarisation with time. Our data lend support to the 
internal shock model whereby the jet velocity increases leading to internal shocks in 
the pre-existing outflow before the jet switches off. The compact nuclear jet is seen to 
re-establish itself within two days, and is visible as core emission at all epochs. The 
energetics of the source are calculated for the possible range of distances; a minimum 
power of 1-10 per cent Lscid is required to launch the jet. 

Key vifords: Accretion, accretion discs - Stars:individual GRS1915+105 - 
Stars:variables - ISM: jets and outflows - Radio continuum:stars - X-rays:stars 



1 INTRODUCTION. 

GRS 19154-105 was the first Galactic source observed t o 
exhibit superluminal motion iMirabel fc Rodri'gued . Il994l) . 
and studies of its radio jets and jet-disc interactions have 
been important in developing our understanding of the link 
between accretion and jet outflows in X-ray binary systems. 
The jets extract mass, energy, and angular momentum from 
the accreting compact object, and in doing so can inject 
significant energy (> 1 per cent of the time-averaged lumi- 
nosity of supernovae) into the ISM (Fender, Maccarone & 
van Kesteren 2005). 

Jet outflows now seem to be ubiquitous in accreting 
systems such as black hole X-ray binaries. Such sources 
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spend the majority of their time in quiescence, with ex- 
tremely low X-ray and radio luminosities and a relatively 
hard X-ray spectrum. In such states, they are thought to 
be 'jet-dominated', such that the power output in the jets 
exceeds that radiated in X-rays (Fender, Gallo & Jonker 
2003). As the X-ray luminosity rises above ~ 7 x 10~^I/Edd, 
the systems become X-ray dominated, with the majority of 
the accretion energy being dissipated as X-rays in the in- 
ner parts of the accretion flow. The jets are then persistent, 
steady, and self-absorbed with a flat radio spectrum (Stir- 
ling et al. 2001; Dhawan, Mirabel & Rodriguez 2000b). At a 
high fraction (~ 10"^) of the Eddington luminosity, the cen- 
tral source makes a transition to a softer X-ray state, passing 
through the so-called Very High/Intermediate states. As this 
happens, the jet velocity increases, leading to internal shocks 
in the steady jet which appear as highly rel ativistic knots 
moving away from the core of the system (e.g. iFender et al.l 
Il999h . At this stage, the steady-state jet outflow is quenched 
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until the system moves back into its low/hard X-ray state 
once more. This unified model is presented for generic black 
hole X-ray binaries by Fender, Belloni & Gallo (20 04), and 
for th e specific case of GRS 1915-^105 by .Fender fc Bellonil 
i2004l) . We describe it in detail, and in the light of the ob- 
servations we present, in §|H1 

The source GRS 1915-1-105 was discovered in 1992 via 
the WATCH instrument aboard the Russian GRANAT mission 
(Castro-Tirado, Brandt & Lund 199 2). The system is be - 
lieved to comprise a K-M III star ijGreiner et all 
and a 14 ± AMq black hole in a 33.5 ± 1.5 d orbit (Greiner, 
Cuby & McCaughrean 2001). The nature of the donor star 
implies that its mass should lie between 1 and 1.5 Mq and 
rules out the stellar wind as the accretion mechanism; the 
mass loss rate from such a star would be insufficient to feed 
the high accretion luminosity of GRS 1915-1-105. Thus it is 
thought that accretion occurs via Roche lobe overfiow. 

In this paper, we present a study of the 2001 March and 
July outburst of GRS 1915-1-105, using data from the Rossi 
X-ray Timing Explorer {RXTE) and the Multi-Element Ra- 
dio Linked Interferometer Network (MERLIN). The X-ray 
observations are discussed in §|5|and the radio observations 
in §121 which are used to place constraints on the jet speed 
and inclination angle, the source distance, and the possible 
deceleration and expansion of the jet knots as they move 
outwards from the core. The radio core emission is analysed 
in §|11 and polarisation maps are presented in We use 
the observations to constrain the energetics of the source in 
§|5| and address the discrepancy between radio observations 
of the jets in this source on different angular scales in §[7| 



2 X-RAY STATES 

The continuum X-ray spectrum of GRS 1915-1-105, in com- 
mon with most other black hole sources, has been modelled 
as a sup erposition of a soft disc b lackbody and a hard power 
law fe.g. lKlein-Wolt et alll2002^ . The blackbody {ksT ~1- 
2keV) is associated with the accretion disc, while the power 
law represents Comptonised emission from a corona, which 
latter component may be ejected leading to a radio flare 
(or oscillations in the case of repeated ejections). The X- 
ray emission is highly variable, which has been interpreted 
as the disappearance and reappearance o f the inner region 
of the accretion disc (iBelloni et al.L[T997al lbh. caused by the 
onset of thermal- viscous instabilities. 



2.1 RXTE Data 

RXTE, launched in 1995, carries two pointed instruments, 
the Proportional Counter Array (PGA) which covers the 
2-60 keV energy range, and the High Energy X-ray Timing 
Experiment (HEXTE) for the 15-250 keV e nergy range. Us- 
ing archival PCA data lBelloni et al] i2000t) found that the 



X-ray variability is the result of transitions between three 
basic states; A, B, and C, categorised according to their 
spectral hardness, temporal variability, and the prominence 
of the thermal disc and power-law components. 

The X-ray data presented here were collected by the All- 
Sky Monitor (ASM) on board RXTE, which scans about 80 
per cent of the sky every orbit. The ASM returns the total 
source intensity in the 2-12 keV band as well as intensities 
in three sub-bands: A (1.3-3.0 keV), B (3.0-5.0 keV) and C 
(5.0-12.2 keV). Two hardness ratios can therefore be calcu- 
lated; HRi ^ B/A and HR2 = C/B. It should be noted 
that the bands and ratios for the PCA differ due to that 
instrument's higher energy range (2-60 keV) and that PCA 
HRi ~ ASM HR2. The ASM count rates in the three bands 
and the HR2 hardness ratio are shown in Figs. and |5| for 
the March and July observations respectively. In all cases, 
the radio outbursts seem to be associated with a rise in the 
X-ray count rate and a degree of spectral softening (a de- 
crease in the hardness ratio HR2). 

PCA observations were made during the March obser- 
vations at the epochs indicated in Fig. and in all cases, 
the temporal variability and spectral hardness showed the 
source to be in State C. PCA observations were also made 
during the J uly outburst, bu t since they have already been 
analysed by IVadawale et alJ i2003l) , we will not reanalyse 
them here. They will be discussed in § 13.21 



3 OBSERVATIONS 

GRS 1915H-105 was observed with the MERLIN interferome- 
ter during two flaring sequences, between 2001 March 24 and 
2001 Aprfl 5, and also between 2001 July 17 and 2001 July 
24. The observing frequency was 4.994 GHz with a band- 
width of 15 MHz. The observing times are listed in Tabled 
The March-April epochs of observation have been labelled 
1-9 sequentially, and the July observations labelled 10- 
16. These target-of-opportunity observations were triggered 
by the ongoing flux monitorin g program at the Ryle Tele- 
scope JPoolev fc Fendeilll997l) . At each epoch, in addition 
to the target source GRS 1915-1-105, observations were made 
of a flux and polarisation angle calibrator, 3C 286, a point 
source calibrator, OQ 208, B 2134-1-004, or B 0552-1-398, and 
a phase reference source, B 1919-1-086, 2.8° away from the 
target source. The point source calibrator used, and the time 
on source for each epoch are also listed in Tabled During 
the observations, the five outstations (Cambridge, Defford, 
Knockin, Darnhall and Tabley) were used, together with the 
Mark 2 antenna at Jodrell Bank. 

The MERLIN d-programs were used to perform ini- 
tial data editing and amplitude calibration, and the data 
were then imported into the National Radio Astronomy 
Observatory's (NRAO) Astronomical Image Processing Sys- 
tem (aips) software package for further data reduction. The 
MERLIN pipeline was then used to image and self-calibrate 
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Figure 1. 15-GHz Ryle Telescope monitoring data, RXTE ASM 
count rates and hardness ratios prior to and during the 2001 
March observations. Vertical dotted lines show the derived zero- 
separation times for the three ejected components, and the top 
set of horizontal lines in the radio and HR2 plots show the times 
of our MERLIN observations. Below those are indicated the times 
of the RXTE PCA observations. Note the spectral softening and 
peak in the count rate at the time of the first ejection event. 
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Figure 2. 15-GHz Ryle Telescope monitoring data, RXTE ASM 
count rates and hardness ratios prior to and during the 2001 
July observations. Vertical dotted lines show the derived zero- 
separation times for the ejected component, and the horizontal 
lines in the radio and HR2 plots show the times of our MERLIN 
observations. 



the phase reference source, and apply the derived corrections 
to the target source, GRS 1915-1-105. The pipeline also calcu- 
lated the instrumental corrections (the D-terms arising from 
signal leakage from right circular polarisation feeds into left, 
and vice versa) using B 1919-f 086, and calibrated the po- 
larisation position angle using 3C 286, assuming a position 
angle of 33° east of north for its electric field vector. Fur- 
ther self-calibration and imaging were then carried out using 
standard procedures within AlPS, using the phase-referenced 
images as initial models for each epoch. The solution inter- 
val was gradually reduced and the number of CLEAN com- 
ponents used in the model was increased until there was no 
further improvement in the images. The images in Stokes 
/, Q and U, were then combined to produce images of to- 
tal intensity, polarisation intensity, and polarisation position 
angle. 

Epochs 1, 6, 7 and 12 were all found to show signifi- 



cant core variability on timescales of tens of minutes over 
the course of the observing run. This made imaging diffi- 
cult, as the sidelobe levels changed with time. For epochs 1 
and 7, the core flux densities were approximately constant 
over a sufficiently long period of time that imaging and self- 
calibration using only a restricted time range was found to 
produce images with acceptably low residual r.m.s. noise lev- 
els. For epochs 6 and 12, the source flux densities were not 
constant for long enough for this approach to work. Instead, 
the data were split into small segments (of length 20 and 30 
minutes respectively; a longer time period was required dur- 
ing the July observations since the source was intrinsically 
fainter). Each segment was self-calibrated to remove any at- 
mospheric jitter, imaged, and the flux density of the core 
fltted and subtracted from the it«;-data. The segments were 
then recombined, and the resulting core-subtracted data set 
was imaged to show any extension. A point source was added 
back in at the fltted core position, with the fltted mean flux 
density of the core, in order to make the flnal images seen 
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in Figs. 1^1^ and [T^ and to enable easier comparisons to be 
made with the other epochs. The core parameters used were 
measured from the original image of the entire dataset for 
that epoch. 



3.1 March outburst 

Fig. El shows a composite of the MERLIN images of the 
March observations. The images have been rotated clockwise 
through 52.5°, so that the southeastern components appear 
on the left, and the northwestern components to the right. 
The core appears to be detected in all epochs (see §2]for fur- 
ther discussion), and we clearly observe components moving 
away to the southeast and the northwest. The southeastern 
components appear to move faster (Fig. |1J and are brighter 
at a given angular separation from the core (Fig. EJ, and 
thus correspond to the approaching jet, while the fainter, 
slower northwestern compone nts are receding from u s. Thi s 
agrees with the findings of 'Mirabej_^^_^x^jrigij^j l )l994l . 
feodn 'gucz & Mirabel (1999) and F ender et alJ l|l99flh . Over 
the course of the 13 days of observation, three distinct south- 
eastern components were seen to be ejected and move out- 
wards, with the middle one being fainter than the first and 
last. These are labelled SE, SE2, SE3 respectively. Owing to 
Doppler deboosting of the receding jet flux density and its 
lower apparent proper motion, only a single northwestern 
component was observed, labelled NW. The fitted flux den- 
sities and angular separations from the core of the different 
components are given in Table|2| 

Assuming that the ejecta move ballistically (i.e. with 
constant velocity), straight-line fits to the angular separa- 
tions can be extrapolated to find the ejection dates of the 
different components. These are given in Table|l]and plotted 
in Fig. 21 and reveal that components SE and NW are con- 
sistent with having been ejected simultaneously. The quoted 
uncertainties take into account both the uncertainties in the 
fitted positions of the core and the jet component, and also 
the uncertainties in the time of observation, taken as half 
the length of th e observing run. The R yle telescope moni- 
toring program dPoolev fc Fendeilll997ll indicates that the 
outburst in the 15-GHz flux density of GRS 1915-1-105 which 
triggered our MERLIN observations peaked at MJD 51990.4 
(see Fig.0. This, and also the observed sharp rise in 15-GHz 
flux density, is in good agreement with our zero-separation 
date of MJD 51989.84 ± 0.42. 
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3.2 July outburst 

IVadawale et alJ l)200,^) carried out a detailed study of the X- 
ray and radio state of GRS 1915-1-105 prior to and during the 
July outburst. The outburst followed a by now familiar pat- 
tern, with a ~ 60 mjy 15-GHz radio flare on MJD 51989.0 
iVadawale et alj . l2003ll preceding a radio-loud plateau state, 
which lasted until the relativistic ejection. They found that 
according to the RXTE PGA observations, prior to July 16, 



Figure 3. Contour maps for the March observations. Solid and 
dashed contours for each image are (v^") and — (\/2") times 
the levels specified on the right-hand axis. The images have been 
rotated clockwise by 52.5° to form the montage. The dotted lines 
correspond to the fitted ejection dates and proper motions of the 
components. The vertical dashed line indicates the core position. 
The beam sizes for each image are plotted in the lower right-hand 
corner. At a distance of dmax = 10.9 kpc, Imas on the imago 
corresponds to a spatial scale of 10.9 au. The component labels 
are shown above the epoch 9 image, and the orientation on the 
sky is also indicated. 
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Table 1. The labels, dates, Modified Julian Dates and on-source times for the 2001 MERLIN observations of GRS 1915+105. 



Epoch 


Date (2001) 


MJD 




Time on source (min) 


Point source calibrator 


\ 


MarrVi 94. 


51992.30225 


it 


2587S 


478.5 




2 




51993.32289 


± 


23779 


438.9 


R 0552+398 


3 




51994.27843 


± 


23574 


433.7 


R 0552+398 


4 


March 27 


51995.29641 


± 


0.25411 


468.0 


B 0552+398 


5 


March 29 


51997.28781 


± 


0.24570 


431.3 


B 2134+004 


6 


March 31 


51999.28593 


± 


0.24383 


447.9 


B 0552+398 


7 


April 2 


52001.27083 


± 


0.25383 


468.0 


B 0552+398 


8 


April 3 


52002.27575 


± 


0.25403 


468.4 


B 0552+398 


9 


April 6 


52005.25733 


± 


0.22043 


408.4 


B 0552+398 


10 


July 17 


52107.98903 




0.23861 


372.1 


B 0552+398 


11 


July 18 


52108.98906 




0.23865 


369.7 


OQ208 


12 


July 19 


52109.98907 


± 


0.23865 


371.9 


B 0552+398 


13 


July 21 


52111.98906 


± 


0.23865 


375.8 


OQ208 


14 


July 22 


52112.98282 




0.23239 


357.7 


OQ208 


15 


July 23 


52113.98278 


± 


0.23236 


361.1 


OQ208 


16 


July 24 


52114.96038 




0.25060 


388.3 


OQ208 



Table 2. The fitted flux densities and angular separations from the core of components for the 2001 March observations of GRS 1915+105. 
The integrated flux densities and positions of the components were found by fitting Gaussians, and the angular separations were all 
measured from the fitted core positions in their respective images. 



Epoch 


MJD 


Core flux 


SE scp. 


SE flux 


SE2 scp. 


SE2 flux 


SE3 scp. 


SE3 flux 


NW scp. 


NW flux 




-51990 


(mJy) 


(mas) 


(mJy) 


(mas) 


(mJy) 


(mas) 


(mJy) 


(mas) 


(mJy) 


1 


2.36 ± 0.20 


38.2 ± 0.4 


51.7 ± 2.2 


31.4 ± 0.04 














2 


3.32 ± 0.24 


24.6 ± 0.3 


75.8 ± 2.4 


21.1 ± 0.03 














3 


4.28 ± 0.24 


29.4 ± 0.3 


98.5 ± 2.5 


19.0 ± 0.3 














4 


5.30 ± 0.25 


27.8 ± 0.5 


112.6 ± 10.1 


10.5 ± 0.6 










52.2 ± 17.0 


2.8 ± 0.4 


5 


7.29 ± 0.25 


19.4 ± 0.2 


152.6 ± 15.6 


1.6 ± 0.2 


50.7 ± 18.3 


2.3 ± 0.3 






53.8 ± 29.5 


1.9 ± 0.3 


6 


9.29 ± 0.24 


49.7 ± 0.2 














96.8 ± 60.3 


4.9 ± 1.2 


7 


11.27 ± 0.10 


26.9 ± 0.4 


234.8 ± 32.1 


2.1 ± 0.5 






50.3 ± 1.9 


51.1 ± 0.4 


151.7 ± 21.1 


1.6 ± 0.4 


8 


12.28 ± 0.25 


19.6 ± 0.3 






192.4 ± 28.0 


2.2 ± 0.4 


65.3 ± 2.4 


26.4 ± 0.3 


115.0 ± 28.9 


1.5 ± 0.3 


9 


15.28 ± 0.25 


94.3 ± 0.3 






240.7 ± 67.7 


0.6 ± 0.3 


155.3 ± 7.7 


4.2 ± 0.3 


167.7 ± 15.6 


1.1 ± 0.2 



the source was in a hard state (the radio-loud C-state, Crl), 
in the xrl class, and then moved to a state of enhanced 
X-ray emission, possibly due to a disturbed accretion disc. 
They were however unable to determine the exact time of 
the X-ray state change. During the radio flare, they saw X- 
ray dips, corresponding to periods of state A, during which 
the Comptonised spectral component of the X-ray emission 
was absent. 

During our July observations, from which the images 
are shown in Fig. |S] the flux densities of the core and the 
jet components were in general much lower than during the 
March observations. The flux densities of the core and SE 
component, and the angular separation between the two, are 
given in Table |H] The only epochs in which an unambigu- 
ous detection of the NW component was made were 12 and 
13. There is marginal evidence in epoch 11 for an elonga- 
tion of the core in the opposite direction to that of the SE 
jet component, which could be interpreted as the receding 
component. There may be a north-western extension to the 
core in epoch 14, though this is not particularly convincing 
(having a flux density of only 1-2 times the r.m.s. noise). 



Again assuming ballistic motion, a straight-line fit to 
the angular separations of the SE component from the core 
gives a proper motion of /iapp = 23.8 ± 2.7 masd~^. The 
fit is good, with x?cd = 0.4, and implies a zero-separation 
date of MJD 52106.1 ± 1.4, which corresponds to July 16 
03:21 UT, in g ood agreement w i th the start of the outburst 
determined by IVadawale et al.l (|200^ . A similar fit to the 
receding component gives ^vec = 11-8 ± 3.5 masd~^, and 
a zero-separation date in agreement with that for the ap- 
proaching component. These two fits are shown in Fig. |S| 



3.3 Proper motions 

For a symmetric ejection event, if the proper motions of cor- 
responding approaching and receding components, /ia and 
respectively, can be measured, it is possible to calculate the 
product 

/3cos6>= (1) 

Ma — A*r 
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Table 3. The fitted flux densities and angular separations of components from the core for the 2001 July observations of GRS 1915+105. 



Epoch MJD-52100 Angular separation Error Core flux density Error SE flux density Error 







(mas) 


(mas) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


10 


7.99 


42.1 


4.1 


28.0 


0.6 


33.6 


0.8 


11 


8.99 


71.1 


6.6 


27.1 


0.4 


8.9 


0.4 


12 


9.99 


90.9 


27.9 


27.7 


0.5 


2.6 


0.5 


13 


11.99 


150.2 


30.6 


9.6 


0.4 


2.0 


0.4 


14 


12.98 


158.1 


70.8 


15.1 


0.3 


0.6 


0.3 


15 


13.98 


179.0 


48.7 


15.1 


0.3 


1.2 


0.4 


16 


14.96 


204.3 


29.7 


16.3 


0.3 


0.8 


0.3 



Table 4. Fitted proper motions and ejection dates of jet components, assuming ballistic motion 



Outburst 


Component 


Epochs 


Proper motion (masd ^) 


Ejection date (MJD) 


^red 


2001 March 


SE 


1-5,7 


21.4 ± 2.0 


51989.84 ±0.42 


1.08 


2001 March 


NW 


4-9 


11.8 ±2.0 


51990.9 ± 2.2 


2.92 


2001 March 


SE2 


5,8,9 


26.8 ±5.9 


51995.4 ±2.6 


0.20 


2001 March 


SE3 


7-9 


27.4 ± 2.5 


51999.6 ± 1.6 


1.47 


2001 March 


SE3 


6-9 


24.7 ± 1.0 


51999.2 ±0.6 


3.09 


2001 July 


SE 


10-16 


23.8 ±2.8 


52106.1 ± 1.4 


0.41 


2001 July 


NW 


0^,12-14 


11.8 ±3.5 


52106.2 ±3.8 


0.11 


1997 Oct 


SE 




23.6 ±0.5'' 




< 1 


1997 Oct 


NW 




10.0 ±0.5'' 




< 1 


1994 Mar 


SE 




17.6 ±0.4= 




< 1 


1994 Mar 


NW 




9.0 ±0.1= 




< 1 



" To better constrain the fit, the zero-separation date ('epoch 0') derived from the flt to the SE component was also used 

*' Data taken from^]endCT^t_al^ ('1999') 

= Data taken from Mirabel & Rodriguc^ ^^9^ 




Observation date (MJD-51980) 

Figure 4. Proper motions for the March observations. The dot- 
ted lines correspond to the best fltting ejection dates and proper 
motions for the different components, accounting for uncertainties 
in both time of observation and measured angular separation. 



where P is the jet speed v/c and 6 is the inclination angle 
of the jet axis to the line of sight. Setting 6 — 90° allows us 
to place a lower limit, /Smin on the intrinsic jet velocity. As- 
suming /3 = 1 places an upper limit, Smax, on the inclination 
angle of the jet axis to the line of sight, and allows an upper 
limit, dmax, to be placed on the source distance, given by 



(2) 



For /ia and fi^ in units of mas d ^ , this may be more conve- 
niently expressed as 

173 



VMaMi- 



kpc. 



(3) 



Assuming that the SE and NW components for the March 
outburst were the approaching and receding components 
from a single symmetric event, then from the fitted proper 
motions given in Tabled we find (3 cos 9 ^ 0.29 ±0.09. This 
gives I3min = 0.29 ± 0.09, 6^^^ = 73.3 ± 5.2°, and d,nax = 
10.9 ± 1.0 kpc. If the distance is specified, then the exact 
values of /3 and 6 can be calculated, as can the bulk Lorentz 
factor r = (1 — P'^)~^^^ , and the Doppler factors of the ap- 
proaching and receding jets, (5app,roc = (r[l =p P cos9])~^ . 
These have been plotted in Fig. |7] 



We can also use the proper motions from the July 
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Figure 5. Contour maps for the July observations. Solid and 
dashed contours for each image are {V2") and -{V2") times 
the levels specified on the right-hand axis. The images have been 
rotated clockwise by 52.5° to form the montage. The dotted lines 
correspond to the fitted ejection dates and proper motions of the 
components. The vertical dashed line indicates the core position. 
The beam sizes for each image are plotted in the lower right-hand 
corner. At a distance of dmax = 10.9 kpc, Imas on the image 
corresponds to a spatial scale of 10.9 au. 



outburst given in Table gl to find PcosO = 0.34 ± 0.14, 
filmax = 70.4° ± 8.5°, and d^ax = 10.3 ± 1.6 kpc. These are 
consistent witii tlie values found for the March outburst, 
so for this reason, and since, owing to the more robust de- 
tection of the receding component, the Marcli values were 
better-determined, the jet parameters for this outburst have 
not been plotted. 

Our derived proper motion s are consiste n t (to within 
errors) with those found bv iFender et al.l ilQQgTl with 
MERLIN and iDhawan et~al\ 1 2000b^ with~the VLB A, but 
grea ter than those found by Mirabel fc Rodriguez! (Il994ll 
and iRodrfguez fc MirabeJ ^1999^ with the VLA. The im- 
plications of this are discussed further in §[7| 




Observation date (MJD-52100) 

Figure 6. Proper motions for the July observations. The dotted 
lines correspond to the best fitting ejection dates and proper mo- 
tions for the different components, accounting for the averaging 
time and the uncertainty in measured angular separation. 



3.4 Source distance 



iMirabel fc Rodrfgued il994ri originally derived a distance of 
12.5 ± 1.5 kpc to GRS1915H-105 from Hi absorption mea- 
surements with the VLA, while their measurements of jet 
knot proper motions constrained the maximum possible 
source distance (corresponding to /3 = 1) to be rfmax = 13. 7± 
0.2 kp c. The proper motion observations of iFender et alJ 
l|l999|) gave a somewhat smaller value, dmax = 11.2±0.8kpc. 
Dhawan, Goss fc Rodriguez (2000a) further constrained the 
distance to GRS 1915+105, again using Hi absorption mea- 
surements. They found d > 6.1 kpc and d < 12.2 kpc, and 
adopted a distance of 12 ± 1 kpc. They also found extra 
absorption along the line of sight to GRS 1915-1-105 com- 
pared with that to the Hll region G 45.46-1-0.06 (at a dis- 
tance of 8.8 kpc), implying d > 9 kpc. Recently however, 
iKaiser et a l.' 1*2004) interpreted the two IRAS sources iden- 
tified by [Rodri guez & Mirabel (.1998,1 as the impact sites 
of the jets of GRS 1915-1-105 on the interstellar medium. 
This would require all three sources to be at approximately 
the same distance from us, 6 < d < 7.4 kpc. They also 
pointe d out that the CO observations of IChapuis fc CorbeJ 
i2004^ allowed the constraint on the distance be i ng be yond 
the G 45.46-1-0.06 to be relaxed. Izdziarski et all feOOSi) dis- 
missed these assertions however, claiming that systematic 
errors on the systemic veloc ity derived from infrared CO 
lines bv lGreiner et al.H2001bll would be small. The systemic 
velocity of l&ein^^t^ together with the Galactic rotation 
curve gave a distance of 12.1 ± 0.8 kpc. Our 3a upper lim- 
its on dmax for the March and July outbursts are 13.9 and 
15.1 kpc respectively. 
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Figure 7. Variation with assumed source distance of jet speeds, 
inclination angles, bulk Lorentz factors and Dopplcr factors for 
the approaching and receding jets. Calculations have been done 
for the best fitting proper motions, the proper motions both low 
by la (minimum values), and the proper motions both high by 
1(T (maximum values). The vertical lines indicate the values of 
dmax, 10.9, 12.5, and 9.5 kpc respectively. 



In summary, the distance would appear to lie between 
6.1 and 12.2 kpc, implying intrinsic jet speeds of 0.57 < P < 
1, from our observations and Fig.|7| But it is not yet possible 
to break the degeneracy between jet speed and inclination 
angle and plac e any more rigo r ous con straints on the intrin- 
sic jet speeds. iDhawan et al.l feOOOlJ) suggested that mea- 
surement of the annual trigonometric parallax (~ 80 mas) 
might in future resolve this discrepancy. Such a measure- 
ment would seem to be of prime importance in deciding 
how relativistic the jets from this source actually are, and 
in constraining the energy budget of the outbursts. 



iFendeil i2003l) remarked that all significantly relativis- 
tic jets lie very close to dmax (which is in fact implicit in the 
definition of dmax since it is the source distance correspond- 
ing to /3 = 1). If the jets in GRS 1915-1-105 are as relativistic 
as has previously b een claimed jMirabel fc Rodrfguezlll994l : 
iFender et aI.L Il99tj l . the source is certainly at a distance of 
more th an 9 kpc. Should it b e as close as 6-7 kpc, as postu- 
lated bv lKaiser et alJ (120041) . the jets, while fast, would not 
be the superluminal ejecta often assumed. 



3.4-1 Are the knots decelerating? 

As the jet knots move outwards, their bulk kinetic energy 
will be thermalised as they interact with their environment. 
Depending on the rate at which this process occurs, it might 
be possible to measure the deceleration of the jet knots, 
as has been done on arcsecond scales for the case of XTE 
J 1550-564 dCorbel et alll2002h . The internal shock model 
of Kaiser, Sunyaev & Spruit (2000) however, predicts that 
for a strong shock (which is believed to be the case), the 
velocity of the shock would be constant. In the case of 
GRS 1915H-105, iMirabel fc Rodn'gued measured sig- 

nificantly slower proper motions for both the approaching 
and receding components than found here, at angular sepa- 
rations of up to 720 mas from the core. However, their data 
were consistent with ballistic motion, so deceleration would 
not seem to explain the discrepancy. 

Fitting our measured angular separations with a 
quadratic function, rather than the straight-line fit required 
by ballistic motion, showed no conclusive evidence for de- 
celeration. The best-fitting values for the three components 
where there were a sufficient number of measurements to 
properly constrain the fits (March SE, NW and July SE) 
were decelerations of 0.7±0.5, 0.4±0.8, and 0.9±1.3 masd~^ 
respectively, all consistent with zero. The x?od values of these 
fits were similar to those for the str aight-line fits, and an F- 
test (e.g. IPfenniger fc Revaj. 1200 5l) showed that in no case 
was the probability greater than 95 per cent, i.e. adding in 
the deceleration parameter was not necessary at any sig- 
nificant level. Furthermore, the quadratic fit to the angular 
separation of the March SE component predicts an ejection 
date of MJD 51990.40 ± 1.5, the date of the 15-GHz fiux 
density peak. Assuming that the core fiare corresponds to 
the ejection of that component, the ejection date would then 
have been prior to the flux density peak, strengthening the 
case against deceleration of that component in our data. 
These data do not therefore support the case for decelera- 
tion. 



3.5 Expansion of jet components 

While no clear expansion of the jet component sizes is seen 
in the images themselves, it should be remembered that the 
images of Figs. |3 and |3 are convolved with the Gaussian 
restoring beam of ~ 75 x 40 mas^. Deconvolving the beam 
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Figure 8. Deconvolved component sizes for components SE and 
SE3 of the March observations, perpendicular (top) and parallel 
(bottom) to the jet axis. Dotted lines show best-fitting linear 
expansions. Labels denoting the epochs of observation are shown 
in the top plot. 



from the fitted sizes to find the intrinsic sizes of the Gaus- 
sian components fitted to the images shows marginal evi- 
dence for the expansion of the jet knots. Fig.|H|shows fitted 
Gaussian knot sizes parallel and perpendicular to the jet 
axis for components SE and SE3 of the March observations. 
Gomponent SE only shows convincing evidence for being re- 
solved perpendicular to the jet axis during epochs 1 and 5, 
whereas component SE3 shows significant evidence for be- 
ing resolved in both epochs 8 and 9. The July observations 
were less conclusive, with the deconvolved component sizes 
initially seeming to increase, but falling off again in the later 
epochs as the fiux density decreased and the knots became 
harder to detect. 

If this component expansion is real, it can be used to 
constrain 0, the opening angle of the jet. We see the mo- 
tion of the jet knots projected on the plane of the sky. This 
means that for a given angular separation, Sxp, the true sepa- 
ration of the knot from the core is ^V/ sin 9. Thus, assuming 
spherical jet knots. 



tan ( 



Csin6' 

A^app 



(4) 



where C, — r/d is the angular size of the source. Without 
knowing the exact source distance, we have only an upper 
limit to the inclination angle of the jet axis to the line of 
sight. Taking only the most reliable point from the evo- 
lution of the SE component, that of epoch 4, gives ^ ~ 
7.4±1.5 masd~^, assuming zero size at the ejection date de- 
rived in §E3 Alternatively, a linear fit to the data (shown as 
the dotted line in Fig.|HJ gives a value ( = 3.7±0.5 mas d~^. 
These two values give upper limits on the jet cone half- 
opening angle of t/imax = 18.3° ± 3.5° and 9.4° ± 1.3° for the 



SE component. Since the receding component corresponding 
to knot SE3 is not detected, 0max is not well-constrained in 
that case. 

At a distance of dmax = 10.9 kpc, an angular expansion 
of 7.4±1.5 masd^^ corresponds to a lateral expansion speed 
of 0.47 ± 0.10c. We note that this is close to the theoretical 
expansion speed of a relativistic plasma, c/\/3 = 0.58c. The 
data are unfortunately not good enough to make a more 
detailed comparison possible in this case however. 



3.6 Flux densities 

Accurately measuring the flux densities of the components is 
a difficult procedure, since in order to adequately sample the 
uv-plsLTie, and hence be able to detect the source structure, a 
long integration is required. But the jet components change 
in both position and flux density with time. During a 12- 
hour observation, the proper motion of 21.4 masd~^ will 
cause a component to move by 10.7 mas. 

For the early epochs where the jet component was suffi- 
ciently bright for its flux density to be determined accurately 
in a short time, the core was subtracted from the data in 
the M?;-plane, and the flux density at the position of the 
extension was measured as a function of time using a di- 
rect Fourier transform of the complex visibilities, binning 
the data into 6-minute time intervals. At flux densities be- 
low ~ 5 mjybeam"^, the noise level in the data made this 
impossible, and the flux density of the extension had to be 
measured by fitting a Gaussian to the extension in the image 
plane. The receding component was at no time sufficiently 
bright to use the direct Fourier transform method. 

Fitting a power-law decay with time to the measured 
flux densities for components SE and SE3, and constrain- 
ing the ejection time to be the extrapolated time of zero- 
separation given in Table |1| gave power-law indices of 

1.8 ± 0.03 and 2.01 ± 0.02 respectively, with Xred values of 

3.9 and 1.7 respectively. The indices are Lorentz invariant, 
so are also applicable to the flux de nsity decay in the knot 
frame llRodrfguez fc Mirabel Il999ll . Since the motion ap- 
pears to be ballistic (§ 13. 4. H . the flux density decays with 
angular separation with the same power-law index. 



3.6.1 Flux density ratios 

The ratio between the flux densities of approaching and re- 
ceding components, Sapp and Srec respectively, measured at 
equal angular separation from the core, is given by 



1 +13 COS d 
l-/3cos6l 



(5) 



where a is the spectral index of the emission, deflned by 
Sv oc I'"", where S,^ is the flux density at frequency v. k = 2 
for a continuous jet, and k = 3 for a jet composed of discrete 
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Figure 9. Decrease in flux density witli angular separation from 
tile core, for tlie Marcii observations. Blaclt points sliow tiie ap- 
proacliing component SE, grey points SE3, and open triangles 
the receding component NW. Dotted black and grey lines show 
the best-fitting power law decays for components SE and SE3. 
Flux densities for the approaching components were measured 
with time, converted to angular separations using the best-fitting 
proper motions for the two approaching components. 



ejecta. However, as pointed out bv lFendeil (|20oJ , a measure- 
ment of the flux density ratio gives no information in addi- 
tion to that already found from the proper motion analysis 
of 13.31 other than the value of the fc-parameter. Neverthe- 
less, we can attempt to estimate this parameter from the 
data. The derived proper motions and zero-separation dates 
of the approaching components were used to convert the 
flux density decay with time found in § 13.61 to flux density 
as a function of angular separation from the core for the 
approaching component. Making a comparison to the flux 
density of the receding component should yield a flux den- 
sity ratio measured at equal angular separation. However, 
Fig. |5| shows that this is not a straightforward procedure 
for these data. At small angular separations, where the flux 
density is high enough to be accurately measured, but the 
angular separations are less accurate, the flux density ratio 
is of order 10-15. At greater angular separations, the errors 
on the flux density determination are high, but the ratio 
is closer to 4-5. It appears that the flux density of the re- 
ceding component is still rising up to an angular separation 
of ~ 100 mas, although the uncertainties on the data are 
large, so the flux density ratio is not constant as a function 
of angular separation. We now turn to a different approach 
to measuring k. 

Miller- Jones, Blundell & DuflFy (2004) derived an ex- 
pression relating the value of /3 cos 9 to the flux density ratio 
in a single image, assuming adiabatic expansion of the jet 
components, with the knot radius scaling linearly with time. 



In such a case, the flux density ratio in a single image is 
given by 



Srec \l~l3cOSeJ 



where p is the index of the electron spectrum, such that 
N{E) — kE^^. We note however that if the rise phase was 
optically thin, and the image used was made at a time when 
the flux density of the receding jet was still increasing while 
that of the approaching jet was decreasing owing to adia- 
batic expansion, then this formalism would be invalid, since 
both knots must be in the adiabatic expansion regime of de- 
creasing flux density. As already remarked. Fig. 1^1 appears to 
show that the flux density of the receding knot is still rising 
up to an angular separation of ~ 100 mas (with a large error 
bar), and only after that does it start decaying, so care must 
be taken when applying the above formalism. 

The only epochs for which corresponding approaching 
and receding jets were both clearly visible and distinguish- 
able from the core were 5, 7, and 13. These gave flux density 
ratios of 0.89 ± 0.17, 1.27 ± 0.47, and 0.76 ± 0.19 respec- 
tively. We note that in two of the three cases, the reced- 
ing jet is brighter than the approaching jet, and the third 
case is in fact consistent with this scenario. This would 
constrain k < p. Since we do not have multiwavelength 
observations, it is not possible to deduce the spectral in- 
dex of th e jet knots, and he nce the value of p = 2q -I- 1. 
However, iFender et al.l (l2002f) observed this same outburst 
at both 4.8 and 8.64 GHz, and found a spectral index 
for the integrated (core plus jet knots) flux density of the 
source which stabilise d at ~ 0.77 ± 0.05 by MJD 51991.9. 
iMirabel fc Rodrfgued 1^9^ measured a spectral index of 
0.84 ± 0.03 for both the approaching and the receding jet 
components (once they had separated from the core) in 
the outburst of 1994 March. Assuming that the spectral in- 
dex was the same in the outbursts we observed, this yields 
k = 2.5 ±0.3 for epoch 5, k = 3.1 ±0.6 for epoch 7, and k = 
2.3 ± 0.4 for epoch 13. These are signiflcantly gr eater than 
the values of 1.3-1.9 found bv lFender et al.l (119991), but con- 
sistent with the value of 2.3 found bv lMirabel fc Rodriguez! 
il994f) . and would seem to imply something intermediate 
between discrete ejections and a steady jet. The high value 
for epoch 7 would seem to argue for a discrete ejection. How- 
ever, it would seem strange for the character of the outburst 
(i.e. the value of k) to change as the knots moved outwards. 
A possible reason for the discrepancy could be that the for- 
mula was not applicable during epochs 5 and 13. The an- 
gular separation of the receding component during epoch 5 
is 53.8 ± 29.5 mas which, from Fig. ^ would imply that k 
cannot be calculated accurately for this epoch using Equa- 
tion 1^1 thus lending more credibility to the value derived for 
epoch 7. Since we have no similar plot for the July outburst, 
it is not possible to evaluate whether or not such an effect 
could also explain the low value of k for epoch 13. 
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Figure 10. Flux density of the core during the 2001 March ob- 
servations. The time range of each observation epoch is indicated 
at the top of the plot. Epochs 1, 6 and 7 all show a fractional vari- 
ation of at least 50 per cent. The black points are for the core, 
and the grey points are for the approaching component, SE for 
epochs 1-5, and SE3 for epochs 7—9. 



4 THE CORE 

Detectable emission seems to be associated with what we 
interpret as the core of the system at all epochs, although 
in some of the observations it appears to be blended with 
an emerging approaching or recedi ng jet component. T his 
is in contrast to the observations of iFender et alj dlQQSTl in 
which the flux density was always dominated by the jet com- 
ponents. Fig. 1101 shows the flux density of the core and jet 
components during the March observations. As will be ex- 
plained further in §|S] the core was not in general found 
to show significant linear polarisation, except possibly for 
epochs 2 and 9 of the March outburst sequence. 

Epochs 1 and 6 show that the core flux density faded to 
very low levels immediately after a flare. This can be iden- 
tifled as the jet s uppression in t he sof t disc-dominated state 
in the model of iFender eraP J2004l) (see § El for details). 
However, by the start of epoch 2, the core flux density had 
recovered and stabilised at a level of approximately 20 mjy. 
We attribute this low-level radio emission to the steady, com- 
pact nuclear jet. This stable level of ~ 20 mJy at 5 GHz was 
also seen following the July outburst. This implies that the 
source has moved back to a harder X-ray state to the right 
of the jet-line (see §|S]for further details) without launching 
a second major ej ection. Only crossin g the line from right 
to left in flg. 7 of IFender et al.1 (l2004ll (hard to soft X-ray 
state) gives rise to the internal shocks and the corresponding 
relativistic ejecta. 

The flare in core flux density seen during epoch 6 of the 
March observations (Fig. IIH corresponds to the ejection of 
component SE3. The core flux density began to rise from 
a flat base level of 9.8 ± 1.2 mJy at MJD 51999.20 ± 0.02, 
peaking at 107.4 ± 1.4 mJy 0.24 d later. Assuming the start 
of the rise phase to be the ejection date of the component, we 
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Figure 11. Flux density of the core during epoch 6 of the 2001 
March observations. This does not follow the fast rise and expo- 
nential decay (FRED) behaviour observed for the integrated flux 
density of such systems. Note the rise from and fall back to very 
low flux densities (< 10 mJy) either side of the core flare. 



get an extra constraint on the proper motion of component 
SE3, used in the flt for epochs 6-9 given in ij l3.2l and Table|l| 

A detailed inspection of Fig. llll seems to show evidence 
for slightly more variable radio emission immediately prior 
to the rise in fl ux density correspon ding to the flare, as orig- 
inally noted bv IFender et alJ (|200J). But we note that this, 
together with the rapid falloff in flux density at the very end 
of the observing run, occurred when the source elevation was 
low (< 15°), at which point the weather begins to affect the 
gain solutions (otherwise good to ~ 5 per cent). A check on 
the flux density of the phase calibrator confirmed that only 
for the time range covering the last two points in Fig. 1111 
were the data thus affected. The rest of the data show that 
the core jet flux density decreased with an e-folding time 
of ~ 0.07 d. The data from epoch 1 show a very low core 
flux density (observed when the source was still at a fairly 
high elevation (~ 35°)). This is 2.5 ± 0.4 d after the de- 
rived ejection date of the SE component, by which time the 
jet knot would have moved out to an angular separation of 
~ 55 mas, i.e. out of the beam. We cannot accurately probe 
the core flux density on its own until this time, as we will 
be unable to decouple the core emission from the jet knot 
emission. But the low flux density seen during epoch 1 sug- 
gests that the steady core jet is not re-established until the 
increase in core flux density 2.6 d after the ejection event. 
Th is is considerab l y longe r than the timescale of 18 h found 
bv lDhawan et alj feoOObI) for the nuclear jet to reform fol- 
lowing the start of a maj or outburst, derived from the date 
of the first observations iFender et alJ . I^999^ of radio oscil- 
lations following the outburst. 

No short-period (20-40 min) o scillations of the core flux 
density such as those described bv lPoolev fc Fenderl (119971) 
were observed. A search was made by measuring the fiux 
density in the uv-plane at the core position in 1-minute 
time bins, having subtracted off any extension from the data. 
The resulting time series was run through one-dimensional 
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CLEAN algorithm (Roberts, Lehar & Dreher 1987), which 
showed no significant peaks in the power spectrum of the 
data on frequencies shorter than 2d^^. 

The flux density of the core was much higher during 
epoch 9, at round 80-100 mJy, and seemed to oscillate with 
a period of approximately 8 hours. We cannot however de- 
termine whether these are true oscillations, since we only 
sample ~ 1.5 of these 'periods' in a 12-hour observing run. 



4.1 Positional accuracy 

The beam sizes for the various images are shown in Figs. |3 
and 13 But the use of phase referencing allows the deter- 
mination of the target position to much greater accuracy, 
set by the uncertainty in the position of the phase calibra- 
tor, the uncertainty in the telescope positions, and atmo- 
spheric phase winding. B 1919-1-086 has a position accurate 
to « 14 mas iBrowne et al.L Il99^ . and being 2.84° from 
GRS 1915-1-105, introduces a positional uncertainty of about 
10 mas if the atmosphere causes phase wraps on timescales 
of about 2 hours (the worst-case scenario in the data). 

The positions of the core component were measured 
prior to self-calibration for each epoch, and a weighted 
mean position was taken for the July observations. This 
gave a B 1950 position of 19'' 12^49" .96941 ± 0=.00009, 
10°51'26".6427 ± 0.0051, with a spread of 0.25ms in RA 
and 15 mas in Dec. Owing to the multiple ejections dur- 
ing the March observations, the core position was harder to 
measure accurately, owing to the blending of the core with 
newly-emerging components. This lead to a much greater 
spread in the measured core co-ordinate, of 2.67 ms in RA 
(1.12 ms neglecting epoch 7) and 17.1 mas in Dec. The 
weighted mean position was 19''12■"49^96945 ± 0^00053, 
10°51'26".6588 ± 0.0060 (B1950). We note that these po- 
sitions all rely on having correctly identified the relatively 
stationary central component in the images as the core. Fur- 
thermore, the systematic errors in the astrometry mentioned 
in the previous paragraph mean that the uncertainty on 
these derived absolute positions is at least of order 15 mas, 
although we did not have an astrometric check source we 
could use to quantify accurately this uncertainty. Since the 
hour angle coverage was approximately the same within the 
each set of observations (March and July), and since the 
same phase calibrator source was used in all observations, 
the errors in the relative positions from epoch to epoch are 
much smaller, and are dominated mainly by the signal-to- 
noise ratio and the atmospheric phase winding. The major 
source of error in the angular separations within a single 
image is the signal-to-noise ratio of the knots. 

The most accurate determination of the position of 
the co re of GRS 1915-1-105 is currently dPhawan et al.L 
l2000b^ 19''15™11."54938±0."00007, 10°56'44."7585±0."001 
(J 2000) on 1998 May 2. The proper motion on the sky was 
determined as 5.8 ± 1.5 mas yr"'^, ascribed to secular paral- 
lax. In the three months separating our sets of observations. 



Table 5. Polarisation parameters for the 2001 observations of 
GRS 1915-1-105. P is the polarised flux density, P/I is the frac- 
tional polarisation, and P.A. is the polarisation position angle. 
Whore there was no detection, the Scr upper limits have been 
given. 



Epoch 


Component 


P (mJy) 




P/I 


P.A. 


1 


Core 


0.67 ± 0.08 


0, 


.017 ±0.002 


54.6 ± 5.1° 


1 


SE 


5.80 ±0.09 


0, 


484 ±0.004 


-67.0 ±0.7° 


2 


SE 


3.59 ±0.10 


0, 


470 ±0.005 


-65.8 ± 1.1° 


3 


SE 


1.63 ±0.09 


0, 


.086 ± 0.005 


-51.9 ±2.3° 


4 


SE 


1.50 ±0.13 


0, 


.143 ± 0.015 


-49.8 ±3.6° 


4 


NW 


0.65 ±0.12 


0, 


.236 ±0.056 


-82.0 ±7.2° 


5 




< 0.32 








6 




< 0.43 








7 


SE3 


5.14 ±0.09 


0, 


.121 ±0.003 


-61.9 ±0.8° 


8 


SE3 


3.54 ±0.10 


0, 


.134 ±0.004 


-59.6 ± 1.2° 


9 


Core 


3.52 ±0.16 


0, 


.037 ± 0.002 


73.1 ± 5.1° 


10 


SE 


2.79 ±0.11 


0, 


.083 ± 0.007 


-90.2 ±3.8° 


11 


SE 


0.63 ±0.09 


0, 


.070 ±0.010 


-77.8 ± 5.8° 


12 


SE 


0.31 ± 0.16 


0, 


.119 ± 0.066 


-74.7 ± 15.1° 


13 




< 0.49 








15 




< 0.49 








16 




< 0.53 









this corresponds to a shift of 1.8±0.5 mas between our March 
and July observations. Our uncertainties are not sufficiently 
small to be able to verify this fig ure. However, comparis on 
between our position and that of iDhawan et all i2000bll is 
difficult, since the two sets of observations used different 
calibrators for the phase referencing, and since MERLIN 
uses B 1950 co-ordinates whereas the VLBA uses J 2000 co- 
ordinates. 



5 POLARISATION 
5.1 Linear polarisation 

Polarisation images were made for all epochs, except epoch 
14 (July 22) , when the polarisation angle calibrator 3C 286 
was not observed. Table|S]details the polarisation properties 
of the detected components. Since the overall flux density of 
the source was relatively lower during the July observations, 
polarisation was not detected at a signiflcant level in the 
SE component after epoch 12, which itself only had a 2a 
detection, and no significant core polarisation was detected 
at any time. The polarisation maps are shown in Fig. ll2l and 
typically have a r.m.s. noise level of ~ 150/iJy beam"'^. The 
SE component appeared to be significantly polarised, at a 
level of between 5 and 20 per cent. In epochs 1 and 9, the 
core was also found to be polarised, and in epoch 4 there 
was marginal evidence for the detection of polarisation in 
the emerging receding (NW) component. 

The polarisation intensity, fractional polarisation, and 
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Figure 12. Contour maps with superposed linear polarisation E 
vectors. Solid and dashed contours are (\/2)" and -(v^)" times 
the 3-(T noise level in the image. The lengths of the polarisation 
vectors are scaled to represent the level of linearly polarised flux 
density P in the image. A vector of length 100 mas corresponds to 
3.33 mjybeam"^. To reduce confusion, only polarisation vectors 
greater than the 3-a level in the total polarisation image have 
been plotted, except in the case of epoch 12, where the vectors 
have been plotted down to the l-a level. The beam sizes for each 
image are plotted in the lower right-hand corner. At a distance 
of 10.9 kpc, Imas on the image corresponds to a spatial scale of 
10.9 au. 
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Figure 13. Linear polarisation parameters for the 2001 March 
outburst. The bottom panel shows the total intensity on the left- 
hand axis (dotted lines) and the fractional polarisation (dashed 
lines) on the right-hand axis. The top panel shows the polarisation 
position angle. The parameters measured for the SE and SE3 
components are shown by squares and diamonds respectively. 



polarisation position angle for the SE and SE3 components 
are plotted in Fig. 1131 As the component moves outwards, 
the polarisation position angle rotates, and both the inte- 
grated and fractional polarisation decrease. If the change in 
position angle is due to changing Faraday rotation as the 
jet component expands and moves outwards, the implied 
changes in rotation measure for the SE and SE3 compo- 
nents of the March observations are 83.3 radm~^ over the 
course of 4 days and 11.1 radm~^ over 2 days respectively, 
whereas the change between epochs 10 and 11 of the July 
observations is 60.1 radm"^ over 2 days. Without multi- 
frequency observations, it is not possible to tell whether the 
observed rotation of the positi on angle is in fact du e to Fara- 
day rotation. But we note that lFender et all (|200^ observed 
a 'rotator' event in 2001 January, when the polarisation po- 
sition angles at 4.8 and 8.4 GHz rotated together, implying 
that Faraday rotation was not responsible. A different mech- 
anism is therefore required at times in this source. Such a 
mechanism could be that the observed jet knot is composed 
of a series of shocks, which fade at different rates, or even 
reshock and brighten. As different shocks come to dominate, 
the observed position angle of the polarisation vector would 
appear to change if the field lines were oriented differently 
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in the different shocks. A more exciting possibility however 
is that we could be seeing genuine rotation of the field lines, 
either tracing an underlying helical pattern in the field l ines 
at large distances from the core fe.g. lGomez et alll200lll . or 
possibly even due to rotation of the knots themselves. But 
with the limited sampling and single frequency monitoring 
available, we cannot distinguish between these possibilities. 



5.2 Circular polarisation 



The ATCA observations of lFender et"ai] l|2002|) found circu- 
lar polarisation of —0.56 ±0.07 mjy (0.32 ±0.05 per cent) at 
4.8 GHz just prior to our observations of epoch 1. Our MER- 
LIN observations put a 3(7 upper limit of 0.5 mJy beam"^ on 
the Stokes V fiux density during epoch 1 when the source 
was at its brightest, although the circular polarisation feeds 
on the MERLIN antennas make Stokes V more susceptible 
to contamination by Stokes / than is the case for the lin- 
early polarised feeds used on ATCA. The ATCA detections 
were interpreted as either intrinsic circular polarisation of 
the synchrotron emission or the conversion of linear to cir- 
cular polarisation in the synchrotron-emitting plasma. In 
either case, a low-energy tail to the electron energy distri- 
bution is implied. 



6 ENERGETICS 

Taking the rise time of the fiare observed during epoch 6, 
At, as a constraint on the volume of the emitting region, 
V = 47r(cAt)'^/3, and assuming that the increase in fiux 
density origi nates from with in this region, then using the 
formalism of iLongaiil 1^9^, we can estimate a minimum 
energy associated with this fiare. 



2/7 / 



VWHz 



4/7 



J, 



(7) 

where V is the source volume, v is the frequency at which the 
luminosity Li, is measured, and {rj— 1) is the ratio of energy 
in protons to that in relativistic electrons. The magnetic field 
corresponding to this minimum energy criterion (close to but 
not identical t o the equiparti tion magnetic field), Bmin, may 
be expressed llLongaiil ^^) as 



1.8 



V 



2/7 



,1/7 



(8) 



In order to evaluate these quantities, the rise time of the 
outburst, the frequency and the luminosity must first be 
calculated in the rest frame of the source. Assuming that 
the fiux density is dominated by the approaching jet, the 
rise time At in the rest frame of the observer is Doppler 
compressed compared to that in the source frame, r, as 



r = 5appAt = — 



At 



/3cos e)' 



(9) 



The rest frame frequency and fiux density may also be found 
for the approaching component, v = 5appJ^', where primed 
quantities denote those in the rest frame of the jet knot, 
and unprimed quantities those in the observer's frame. Since 



is a Lorentz invariant, 



Si 



c3+a 
'-'app 



S, (r[l-/3cos6']) 



3-l-Qi 



(10) 



where the extra power of a corrects for a source spectrum 
not being fiat. The monochromatic luminosity is expressed 
as = 47rd'^5'^. Since this depends on the distance to the 
source, the inferred minimum energy and minimum energy 
magnetic field are plotted as a function of source distance in 
Fig. 1141 These are derived under the assumptions that the 
filling factor of the source is ~ 1, that the spectral index 
of the emission is a = 0.75 (consistent with the spectral 
index for the i ntegra ted emission derived for this outburst by 
iFender et al.l (1200211 ') . that the radio emission extends over 

a range in frequency such that i^mix"^'''^ <^ '^mi'n~^'''^i ^.nd 
that the observing frequency f = Vmin, the lowest frequency 
down to which radio emission is seen. This puts a lower limit 
on the minimum energy of the outburst. The low-energy tail 
to the electron energy distribution implied by the circular 
polarisation observations outlined in § 15.21 would however 
significantly increase the total energy requirement. 

The minimum energy is a more slowly- varying function 
of source distance than the minimum power, 2i?min/5app At 
(the factor of two arising from the assumption of symmetric 
approaching and receding components), since the Doppler 
factor 5app decreases rapidly as the source distance ap- 
proaches dmax (see Fig.0. 

Knowing the luminosity of the source (in its rest frame) 
and the minimum energy field, we can estimate the total 
number of relativistic electrons. For a power-law electron 
spectrum, N{E) = kE~^ dE, assuming that the maximum 
electron Lorentz factor is large, we can express 

K = (p - l)ntot(meC^)*'"SmTn . (H) 

where ntot is the electron number density, rUa is the electron 
mass, and 7min the minimum electron Lorentz factor. We 
assume that each electron radiates at frequency 



(12) 



where e is the charge of an electron. Assuming that the min- 
imum Lorentz factor electrons correspond to the o bserving 
freque ncy, we can then use equation (19.17) from iLongaiJ 
ljl994h to derive 



iVtot = 



B 



A(a){p — l){mcc 



2^p-l 



27rmc 



(13) 



where A'tot i s the total num ber of relativistic electrons. A{a) 
IS given by (|Longaiii Il99l p. 292), and is equal to 594 for 
Q = 0.75. 



As explained by iFender et al.l (^9^, the kinetic en- 
ergy associated with the bulk motion of the jet knot is 
given by (F — 1)-E'min for an electron-positron jet, or by 
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Figure 14. Variation with assumed source distance of minimum 
energy {-Eniin)i minimum energy magnetic field (-Binin)i minimum 
power requirement (Pmin)i total number of relativistic electrons 
(A'tot), and minimum kinetic energy of jet knot {K^i^). Calcu- 
lations have been done for the best fitting proper motions, the 
proper motions both low by la (minimum values), and the proper 
motions both high by Icr (maximum values). The vertical lines in- 
dicate the corresponding values of dmax for these cases, 10.9, 12.5, 
and 9.5 kpc respectively. The dotted line in the minimum power 
plot shows the Eddington luminosity for a IAMq black hole. 



(r — l)(i5min + A'totJTipC ) for SLU electron- proton jet with 
'cold' protons, where mp is the proton mass. These two cases 
are plotted in the top panel of Fig. 1141 and, unless d is very 
close to dmax, differ only by a factor 2-3. 

The uncertainties in the distance set out in i) l3.4l should 
be borne in mind when considering the values in Fig. 1141 
The minimum energy, minimum power, and kinetic energy 
are all lower limits, which would rise if the source was far 
from equipartition, which is certainly possible given that the 
jet knot is decaying and expanding, and is thus manifestly 
not in a steady state. Furthermore, all these estimates would 
be further increased if, as seems to be the case from the de- 
tection of circular polarisation, the electron energy distribu- 
tion carries on as a power-law down to energies significantly 
lower than those associated with the observing frequency via 
Equation ^] 



7 PROPER MOTION DISCREPANCY 

These observations seem to have confirmed the observed 
discrepancy between the measured proper motions on mil- 
liarcsecond scales, as measured by the VLBA and MERLIN, 
and those measured on arcsecond scales with the VLA. The 
proper motio ns for the SE compone nt of the 1997 October 
flare found bv lDhawan et alj ^2000b^ over the course of 2.5 h 
correspond to 22.1 ± 1.9 masd~^, consistent with the value 
of 23.6=b 0.5 masd"^ found with MERLIN for the same event 
iFondor et al.L |l 999|). and im p ly no d eceleration between 50 
and 320 mas. iDhawan et al.l (l2000b^ also imaged discrete 
ejecta for an event in 1998 May, finding a proper motion of 
22.3 ± 1.7 mas d~^ for th e approaching compon e nt ov er the 
course of 4.5 h. However. iRodrfguez fc MirabeJ (Il999f) con- 
sistently measured proper motions of ~ 17 masd~^ with the 
VLA for four different o utbursts in 1994 January-Ap ril (in- 
cluding that analysed bv lMirabel fc Rodn'guezl . ll994l) . They 
saw no evidence for deceleration between 80 and 1420 mas. 
In this paper, we again measure higher proper motions with 
MERLIN on scales of 50-300 mas, with no evidence for de- 
celeration. 

The discrepancy between t he VLA and MERLI N mea- 
surements was attributed by iFender et alJ il999ft either 
to intrinsic differences in the jet speeds for the two out- 
bursts, or to resolution effec t s bet ween the two arrays, 
whereas IRodrfguez fc Mirabell (Il999ll also suggested a pos- 
sible change in the angle of ejection. Our observations sug- 
gest that the jet velocity does not vary dramatically between 
outbursts, since the proper motions measured for both the 
2001 March and July flares were consistent with one an- 
other, and with those measured bv lFender et alJ il999l) and 
[bhawan ct al. (200!^. However, we note that with such 
a small sample size, this conclusion cannot be regarded as 
definitive. 



As well as the differences in angular resolution between 
these two seemingly discrepant sets of results, there is also a 
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difference in the observation dates. All the VLA observations 
were taken in in 1994, while the MERLIN and VLBA data 
were taken in 1997, 1998 and 2001. This raises the possibil- 
ity that the discrepancy could be a result of jet precession. 
If the jet angle had precessed towards us between the 1994 
and the 1997 observations, the measured proper motions of 
the approaching component would be greater owing to rel- 
ativistic effects. However, those of the receding co mponent 
would then be correspondingly lower (see fig. 4 of iFendeil 
l2003l) . in contradiction to the observations. The VLA mea- 
sured receding proper motions of 7-9 mas d~^, whereas the 
MERLIN observations give 10-12 mas d~^. Also, we might 
expect to observe a similar change between the 1997 and 
2001 observations as seen between 1994 and 1997, but the 
former two results are in fact consistent with one another. 
These arguments would seem to rule out the precessio n the - 
ory, and lend weight to the suggestion ofj^ndcr ct aL| (^g^) 
whereby the lower angular resolution of the VLA resulted in 
the blending of components and hence lower proper motions. 
If this is the true explanation however, it begs the question 
of whether MERLIN itself is resolving the components ade- 
quately, and whether with further improved resolution, one 
would measure still faster proper motions. Certainly none of 
the MERLIN observations show closely-spaced components 
that might be blended at the VLA resolution of 200 mas. 
That all the VLA observations should consist of such super- 
positions of ejecta, whereas none of the MERLIN observa- 
tions did seems to be too much of a coincidence. Until this 
issue is resolved, the derived jet speeds should be consid- 
ered open to possible future revision upwards, making this 
source even more energetic than is currently believed. Ob- 
servations of further outbursts with multiple arrays seem to 
be required to properly resolve this issue. 

A different explanation could be provided by the en- 
vironment of GRS 1915-1-105, and the medium into which 
the jet knots propagat ed. The source was o r iginal ly discov- 
ered in 1992 August iCastro-Tirado et alj . Il992f) and has 
not since returned to full quiescence. B ATSE was launched 
in 1991 April, and lPaciesas et al] il996l) state that no flares 
were detected prior to 1992 May 6. If prior to its discovery 
in 1992 the source had been in its true quiescent state with 
a very weak jet, the environment of the source would have 
been the relatively dense, low-velocity interstellar medium 
(ISM). An outburst at this time would have quickly run into 
the ISM, and the interaction between the jet material and 
the ISM would have determined the advance speed of the 
external shocks thus generated. We note that the infrared 
jets in GRS 1915-1-105 detected by Sams, Eckart & Sunyaev 
(1996), which have never since been detected, could corre- 
spond to the heating of the surrounding ISM by the jet. 
As the jet/ISM boundary moved further out over time, the 
jet would have been less energetic at the boundary, and 
unable to heat the surrounding ISM to detectable levels. 
Prior to later outbursts, the plateau-state steady jet could 
have inflated a bubble in the surrounding ISM similar to 
that seen in Cygnus X-1 (Gallo et al., in prep.), evacuating 
the region through which the jet propagates. Consequently, 
in later outbursts we see only the emission from the inter- 



nal shocks as the high Lorentz factor jet catches up with 
the lower- velocity steady jet, since the jet/ISM boundary is 
considerably more distant from the core of the system. Fu- 
ture VLA observations should therefore measure the same 
proper motions as those measured by MERLIN and the 
V LBA. This scenario would imply that the lobes identifled 
bv lKaiser et al.l (|20o3) are unlikely to be physically associ- 
ated with GRS 1915-1-105, strengthening the case for it to 
lie at a distance of ~ 11 kpc. We reiterate the need for fur- 
ther VLA observations of the relativistic ejecta during an 
outburst of GRS 1915-^105. 



8 JET EJECTION 

IVadawale et alJ i2003f) suggested that the observed flares of 
GRS 1915-1-105 with discrete 'superluminal' ejecta occurred 
after a radio-loud hard X-ray state (the Crl state) in which 
a continuous flat-spectrum jet existed in the radio pre-flare 
plateau state. They postulated that an X-ray dip into state 
A corresponded to the ejection of the Comptonising cloud in 
the core of the system at higher velocity than that of the pre- 
existing continuous jet, creating an internal shock in the flow 
which was observed as discrete, fast-moving, radio-emitting 
ejecta. This model is supported by observational evidence, 
since Gallo, Fender & Pooley (2003) had previously noted 
that steady jet outflows for generic black hole X-ray binaries 
were signiflcantly less relativistic than transient outbursts, 
whereas in the speciflc case of GRS 1915-1-105, lower lim- 
its on the sp eed of the compact je t in the plateau state 
were found bv lDhawan et al ] ^2000b^ and Ribo, Dhawan & 
Mirabel (2004), who measured speeds of /3 = 0.1 and 0.3-0.4 
respectively. 

iFender et al.l il997f) made the connection between in- 
frared and radio flares, suggesting that the two had 
a common origin in synchroton emission from a jet. 
lEikenberrv et all 1^9^ studied the simultaneous infrared 
and X-ray behaviour of the source, and found a close link 
between X-ray and infrared flares, including evidence that 
an X-ray prec ursor spike was assoc iated with the start of 
a flare. Finally, iMirabel et al.l il998t) made simultaneous X- 
ray, infrared and radio observations of GRS 1915-f 105 at two 
epochs in 1997, and found that the infrared/radio flares were 
associated with hard X-ray dips, also suggesting that the 
ejection event occurred at the time of an X-ray spike at the 
end of the spectral softening. 

iFender fc Belloiiil l|200^ created a synthesis of the 
above pictures, together with other work, to suggest that 
above a certain X-ray hardness, the source produces a steady 
outflow. As the X-ray spectrum softens, the jet velocity in- 
creases monotonically until it crosses some 'jet-line', where 
the soft X-ray flux peaks. At this point, the jet velocity 
increases rapidly before the jet is shut off, and the faster- 
moving material causes internal shocks as it collides with 
the pre-existing flow. These shocks then appear as discrete 
relativistic ejecta with high bulk Lorentz factors, such as 
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those we have observed. GRS 1915+105 appears to be con- 
tinually crossing and recrossing this jet-li ne, producinK the 
oscilla tion events such as those reported bv lPoolev fc Fendeil 
il997l) as the jet switches on and ofT. If the source spends 
sufficient time in the X-ray hard state C, there is enough pre- 
existing material in the continuous jet for internal shocks to 
give rise to the discrete, highly relativistic ejecta seen in ra- 
dio observations. This scenario is discussed in more detail 
bv lFender et alj i200'j) . who discussed the energetics of the 
internal shock model, finding that 5-40 per cent of the ki- 
netic energy of the fast jet could be released in the shock, 
which increases the energy budget of §|S]by up to an order 
of magnitude. A physical explanation of the above scenario 
was also proposed, whereby the inner edge of the accretion 
disc moves closer to the compact object as the accretion rate 
rises. This is supported by the increase in the fre quency of 
the lo w-frequency QPOs (figs. 5, 8, and 9 of Vadaw ale et al.L 
|20o3), and the values of the inner disc radius derived from 
model-fitting to X-ray spectra. The escape velocity from this 
point therefore rises, increasing the jet Lorentz factor until 
an internal shock is formed as the jet switches off, the reason 
for which is not yet entirely clear. Observational evidence for 
the jet switching off can be seen in Fig. 1101 where the core 
flux density falls to very low levels in epochs 1 and 6 in the 
immediate aftermath of the ejections of components SE and 
SE3 of the March outburst sequence. 

We see from Fig. ll4l that the minimum power involved in 
launching the jet is of order 1-10 per cent of the Eddington 
luminosity for a 14M0 black hole, unless the true distance 
to the sourc e is very close to th e calculated dmax, when it is 
even greater. lOallo et alJ l|200,'l l noted that in low/hard state 
black holes, jet suppression occurs at a constant fraction of 
a few per cent of the Eddington rate as a source moves into 
a high/soft state. Thus it appears that the jet could be the 
dominant power output channel during the last moments 
before it switches off. 



shocks propagate through the steady jet unaffected by the 
external medium. 

We have observed linear polarisation arising from the 
approaching jet component which decreases with increasing 
angular separation from the core, with a very gradual ro- 
tation of polarisation position angle with time. With the 
current data, we are unable to distinguish between the pos- 
sibilities of Faraday rotation, brightening and fading of dif- 
ferent shocks within the beam, rotation of the knots, or a 
genuine underlying helical structure in the magnetic field. 

We have demonstrated that the energetics of the sys- 
tem, often calculated for an assumed source distance, vary 
substantially with source distance, which could lie anywhere 
between 6 and 12.2 kpc according to current estimates. How- 
ever, the minimum power required for the outburst is 1-10 
per cent of the Eddington luminosity for the 14Mq black 
hole believed to lie at the centre of the system, although 
this may be more if either the source lies very close to dmax 
or does not satisfy the minimum energy criterion. 

The data provi de support f or th e internal shock 
model proposed bv l Kai ser_et_ai] <|2QQ0|) and refined by 
IVadawale et alJ (120031 and iFender et all ^0^), whereby 
the jet velocity increases rapidly at the start of an outburst, 
before shutting off after the Lorentz factor of the ejected ma- 
terial peaks. The increasing velocity causes internal shocks 
which light up the underlying outflow and appear as discrete 
ejecta moving outwards from the core with constant veloc- 
ity. Shocks are only produced once, on crossing the 'jet line' 
from a hard to a soft X-ray state; after the initial ejection, 
no second set of shocks is seen prior to the source moving 
back to the X-ray state C with steady nuclear jet emission. 
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9 CONCLUSIONS 

We have observed two flaring sequences from GRS 1015+105 
in 2001 March and July. The March sequence showed multi- 
ple ejection events, whereas only a single pair of ejecta was 
observed in July. We measured proper motions of 21.4 + 2.0, 
24.7 ± 1.0 and 23.8 + 2.8 masd"^ for the three bright ap- 
proaching components, and 11.8 + 2.0 and 11.8+3.5 masd~^ 
for the two detected receding components, consis tent with 
those found with MERLIN bv lFender et alJ The 3a 

upper limits on the source distance are 13.9 kpc for the 
March outburst and 15.1 kpc for the July event. 

We propose a possible explanation of the discrep- 
ancy between the proper motions measured with the VLA, 
whereby prior to 1992, the source had been in quiescence, 
such that the environment surrounding it had been filled in. 
Subsequently, the jet has inflated a cavity so that internal 
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